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The notion of fluctuation-induced forces is generalized to the cases where the fluctuations have 
nonequilibrium origin. It is shown that a net force is exerted on a single flat plate that restricts 
scale-free fluctuations of a scalar field in a temperature gradient. This force tends to push the object 
to the colder regions, which is a manifestation of thermophoresis or the Soret effect. In the classic 
two-plate geometry, it is shown that the Casimir forces exerted on the two plates differ from each 
other, and thus the Newton's third law is violated. 
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I. INTRODUCTION 



In his pioneering work in 1948, Casimir introduced the 
notion of fluctuation-induced forces when he showed that 
quantum fluctuations of the electromagnetic fields can 
generate measurable long-ranged forces between conduct- 
ing plates [1]. The idea, however, has been subsequently 
generalized to various cases where: (i) the fluctuations 
are of classical thermal origin, (ii) the fluctuating media 
are complex fluids, and (iii) the boundaries are moving 
[2-5]. _ 

In light of the extensive development of the original 
idea of fluctuation-induced interactions, it may seem 
quite natural to ask what happens if the fluctuations that 
cause the interaction are of nonequilibrium nature. A 
first attempt towards answering this question was made 
by G.I. Taylor already in 1951, when he was studying the 
swimming mechanism of flagella (tails of spermatozoa) 
and the theoretically challenging question of the possibil- 
ity of "swimming at low Reynolds number" [6]. In this 
classic work, Taylor has shown that two flagella under- 
going wavelike harmonic deformations have a tendency 
to attract each other, because their undulations intro- 
duce (nonequilibrium) hydrodynamic fluctuations in the 
surrounding viscous fluid that could induce an effective 
interaction between them. More recently, another gen- 
eralization has been introduced where localized sources 
drive the fluctuating medium out of equilibrium [7]. We 
also note that the idea of depletion forces [8] , which are 
also fluctuation-induced in origin, has been generalized 
to systems driven out of equilibrium [9] . 

A particularly interesting way of driving a thermally 
fluctuating system out of equilibrium is by imposing spa- 
tial or temporal variations in the temperature profile. 
Many interesting phenomena are known to appear under 
such conditions, ranging from anomaly in diffusion to 
the old problem of thermophoresis or Soret effect [10,11]. 
In this effect, whose microscopic mechanism is still sub- 
ject to debate [12], solutions of colloidal particles [13] or 
polymers [14] that are placed in a temperature gradient 
experience separation, which signals the appearance of 
a net driving force on particles as a result of the tem- 



perature gradient. It has also been recently shown that 
inhomogeneity in the temperature profile could lead to 
deviation of the self-diffusion coefficient of a tracer un- 
dergoing Brownian motion from the simplistic local Ein- 
stein's relation [15]. 
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FIG. 1. (a) A single plate immersed in a fluctuating 
medium, (b) Two parallel plates separated by a distance H 
in a same medium. In both cases II represent the stresses 
exerted on a plate from the right- or the left-hand side of it. 

Here, we attempt to generalize the notion of Casimir 
forces for nonequilibrium systems in which tempera- 
ture profile is not uniform [16]. We consider a near- 
equilibrium fluctuating system in which we can define 
local and instantaneous temperature. We study the dy- 
namics of a classical scalar field governed by a Langevin 
equation, in which the strength of the noise is propor- 
tional to the local and instantaneous temperature in the 
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system. We consider external objects that restrict the 
fluctuations of the scalar field and calculate the resulting 
fluctuation-induced stress on these objects. We concen- 
trate on the examples of a single flat plate as well as two 
parallel plates, as depicted in Fig. 1, immersed in vari- 
ous temperature profiles. For the case of a single plate, 
we find that a net force is exerted on the plate due to 
the imbalance of the fluctuation-induced forces from the 
two sides, which has a tendency to push it to the colder 
regions. We also calculate the Casimir force between two 
parallel plates in nonuniform temperature profiles and 
find that the two plates experience different forces in vi- 
olation of the third law of Newton. 

The rest of the paper is organized as follows: In Sees. 
II and III, we introduce the Langevin dynamics of the 
scalar field, which we expect to capture the most phys- 
ically relevant consequences of nonequilibrium tempera- 
ture profiles, and define the quantities of interest that are 
to be calculated. Section IV is devoted to the discussion 
of the example of a single plate placed in various tem- 
perature profiles, which is followed by the corresponding 
discussions for the two-plate geometry in Sec. V. Finally 
Sec. VI concludes the paper, while an Appendix outlines 
some details of the calculations. 



II. THE MODEL 

Let us consider an equilibrium system described by a 
scalar field 0(R, t) that undergoes thermal fluctuations 
in a space that is bounded by a number of external ob- 
jects. The field could represent a component of the elec- 
tromagnetic field (e.g. the electric potential) in a dielec- 
tric medium or vacuum [1-3], an order parameter field 
for a critical binary mixture or a magnetic system [4], a 
massless Goldstone mode arising from a continuous sym- 
metry breaking such as nematic liquid crystals or super- 
fluid helium, an clastic deformation field for fluctuating 
membranes and surfaces, or the electrostatic potential in 
charged fluids at very low salt concentrations [5] . In all of 
the above systems, it is possible to write down a Hamil- 
tonian for the fluctuations around the equilibrium state, 
which in Gaussian approximation reads 

H = y J d 3 R [V<HR)] 2 , (1) 

where K is an elastic modulus describing the stiffness of 
the system for the fluctuations around equilibrium state. 

We consider the situation where temperature has a 
slowly varying profile such that thermal equilibrium can 
be achieved locally. This means that we have many heat 
reservoirs locally in contact with our system, which have 
different temperatures. The temperature difference be- 
tween neighboring reservoirs should be sufficiently small 
so that the condition of local equilibrium is fulfilled. Sim- 
ilarly, we can consider a temperature profile with tempo- 
ral variations that are sufficiently slow such that instan- 
taneous equilibrium state can be defined. We would like 



to drive the system described by the Hamiltonian in Eq. 
(1) out of equilibrium by imposing such temperature pro- 
files and examine the corresponding fluctuation-induced 
forces exerted on external boundaries in the medium. 

Our specific choices of boundaries are sketched in 
Fig. 1. We consider both the case of a single plate of 
area A, as well as two such parallel plates separated by 
a distance H. We choose the z axis to be the normal 
direction to the plates, such that the three dimensional 
position vector can be represented as R = (r, z) . We 
assume the Dirichlet boundary condition, which means 
that the fluctuating field 4> is restricted to vanish on the 
plates. 

The physical quantity that we would like to calcu- 
late for such systems is the pressure or the normal force 
per unit area exerted on each plate from the fluctuating 
medium. In the prescribed geometry, this will be U zz — 
the zz component of the stress tensor. Using the Hamil- 
tonian in Eq. (1), one can calculate the components of 
the stress tensor in terms of the two point correlation 
functions of the field. For example, the local pressure 
exerted from the right side on the plate that is located 
at z — is given as 

n+ (r, t) = ~d z d z , (0(r, z, t)<j>{r, z', *)) \ z . z ,^ 0+ , (2) 

where the minus sign denotes that this pressure tends to 
push the plate to the — z direction. (Note that the nor- 
mal to plate at this point is in +z direction.) A similar 
definition and expression can be given for the pressure 
from the left side II ~ z , and the total normal force per 
unit area 

n tot = n--n+, (3) 

exerted on the plate in the positive z direction. Note that 
due to the translational symmetry of the infinite plates 
in the parallel directions, the lateral forces are all zero, 
even for the case where temperature is not uniform. 

III. DYNAMICAL RELAXATION TO LOCAL 
EQUILIBRIUM 

For calculating the correlation functions in the local 
and instantaneous near-equilibrium states resulted from 
a prescribed temperature profile, we use the fact that the 
such states could be reachable through an equilibrium 
dynamical relaxation. We assume a dynamical Langevin 
equation as 

7d^(R, *) = KV 2 <P(R, t) + r/(R, t), (4) 

with the boundary conditions <p(r, 0, t) = <f>{r,H,t) = 0. 
Here 7 is the friction coefficient of the system, and the 
random force ry(r, t) is a Langevin noise whose spectrum 
is given by the local and instantaneous fluctuation — 
dissipation theorem: 
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(r?(R, t)iy(R', i')) = 2 7 /c B T(R, i)<5 3 (R - R')<J(t - i 7 ), (5) 

where T(R, t) is the temperature profile in the system. 

For solving the above Langevin equation, we define the 
appropriate Green's function as the solution of the fol- 
lowing equation 

{■yd t - if V 2 )G(R, t; R', t') = <5 3 (R - K')S(t - t'), (6) 

with the boundary conditions 

G(R, t, R', f )|r on boundary = 0, (7) 

for each value of R'. The solution of the Langevin equa- 
tion can then be written as 

0(R, t) = j ^R'dt' G(R, t; R, i') ^(R, i'), (8) 

which yields the two point correlation function of the field 

as 

(<£(R,t)<£(R / ,t / )) =2 7 fc B y d 3 Ridti T(Ri,ti) 

x G(R, i; Ri, ti) G(R, i'; Ri, (9) 

Note that the fluctuations in the 0-field are affected by 
the entire temperature profile through the scale-free re- 
laxation dynamics governed by Eq. (4). 

To obtain the form of the Green's function for the spe- 
cific geometries discussed above, we proceed by direct 
solution of the differential equation. We focus only on 
the two-plate case, since the single-plate Green's func- 
tion can be deduced from it by simply setting H = 0. 
Because of the translational symmetry in time, as well 
as the space direction parallel to the plates, we introduce 
the Fourier transformation 

G q , w (z, z) = J d 2 vdt (&■('-* )-*•>«-*') G ( rj Zj t . r ' 7 z ' ; 

(10) 



which satisfies the following equation: 

{Q 2 _-d 2 z )G^{z,z') = ^5{z-z l ) 



(11) 



where Q± = \/q 2 ± i-fiv/K. We can devide the space 
into three distinct subspaces: left (L), middle (M), and 
right (R), as shown in Fig. 1(b), and label the Green's 
function in each subspace with the corresponding indices: 
L, M, or R. The boundary conditions can be summarized 
as 

G^(0,z') = G^ R (L,z') = 0, 
G^J-oo,z') = G R (oo,z') = 0, 
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G^ R (z>-,z') = G^ R (z'+,z>), 



d z G^' R (z'-,z') - d z G^' R (z'+,z')) = 1. (12) 



Using the above boundary conditions we can solve 
Eq.(ll) and directly obtain the closed form for Green's 
function in the different parts of the space as 



Gq,Lo( Z > Z ') 



G^(z,z') 



G R ^(z, z') 



-sinh[Q_z] 2M|_£l ; z>z >, 
-sinh[Q_z'] ; z < z' , 



sinh[Q_z'] sinh[Q- (z-H)] 
KQ- sinh[Q_i?J ' 

sinh[Q_z] sinh[Q_ (z'-H)] 
KQ- sinh[Q_i?] ' 

sinh[Q_(z' - H)] ° xph g ( r g)1 ; * > 



z > z' 



z < z' 



sinh[Q_(z-F)] 



cxp[-Q_( z '-g)] 
-R"Q_ 



Z < Z ', 

(13) 



The above expressions for the Green's function can 
be used to calculate the normal stress exerted on the 
plates for arbitrarily given temperature profiles. While 
the equilibrium fluctuation-induced forces are universal, 
it appears that forces induced by nonequilibrium fluctu- 
ations are not universal and do depend on the specific 
choice of dynamics as well as the microscopic details, as 
will be shown in the next two sections. 



IV. SINGLE-PLATE GEOMETRY 

In this section we focus on the example of a single 
plate that is embedded in a fluctuating background with 
local temperature. When the temperature profile is uni- 
form, the stress resulted from the field configurations on 
the left side cancels completely with the corresponding 
counterpart coming from the right side, and the plate 
experiences no net force. However, when temperature 
is not uniform the strength of the fluctuations are also 
not uniform, and we expect a force imbalance due to the 
asymmetry caused by the temperature gradient. To fur- 
ther simplify the problem, we assume that temperature 
only depends on the z coordinate, i.e. it is independent 
of time and the parallel coordinates. The two contribu- 
tions to the normal stress at the position of the plate can 
be written as 



n --Wf^ />-<=> ^ +a -'=- 



(14) 



and 

n+ = 



7 fe B 

" K 



J WF k 



dzT(z)e-( Q + +Q ->. (15) 



Due to the nonuniformity of the temperature profile, the 
nonuniversal contributions to the stress on the two sides 
do not cancel out. Therefore, the forces mediated by 
nonequilibrium fluctuations are not universal and they 
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also depend on the specific imposed temperature profile. 
Examining the behavior of Eqs. (f 4) and (15) above, we 
can distinguish between three different categories: (i) the 
case where temperature is continuous on the plate and 
its asymptotic values are equal at z = — oo and z = +00, 
(ii) the case where temperature is continuous on the plate 
but its asymptotic values at the two infinities are not the 
same, and (iii) the case where the temperature is not 
continuous on the plate. Each category is examined sep- 
arately below. 



n to t = c\ 



k B (T L -T R ) 



(17) 



where a is a microscopic length scale, and c\ is a nonuni- 
versal numerical constant of order unity. The above re- 
sult can be rewritten as 



n. 



VT, 



(18) 



where VT = Tn 2 ^ L > leading to a Soret coefficient [17] 
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(19) 



Note that the length scale £ appears in an anomalous 
way, in that it is not represented entirely in the form of 
VT. 





FIG. 2. Cross section of temperature profile for different 
categories, (a) Temperature is continuous on the plate but 
its limiting values at ±00 are not equal, (b) Temperature is 
continuous on the plate and its values at ±00 are equal, (c) 
and (d) Temperature is not continuous on the plate. 



B. Continuous on the plate and T(— 00) = T(oo) 

As another example, which represents the case when 
the temperature is continuous on the plate and has equal 
asymptotic values on the far left and right sides of the 
plate, we choose [see Fig. 2(b)] 



T{z) 



T - (5T (e~ 2 /« - e"" l2 /«) 
T + 5T (e z /Z - e mz /t) 



z > 0, 
z < 0, 



(20) 



for any arbitrary value of to. The above temperature 
profile yields 



n to t = c 2 (m - 1 



7- 9 ' 



(21) 



where C2 is a numerical prefactor of order one. The above 
result can be recast in the form 



ntot Ivt, 



(22) 



A. Continuous on plate and T(— 00) 7^ T(oo) 

We assume a specific profile for the temperature as [see 
Fig. 2(a)] 



Tl ; z < 

T{z) = { -\{T L - Tr)| + \{T L + T R ) ■ -e<^<C, 
Tr ;i< z, 

(16) 

which exemplifies the general class where the two half- 
spaces are kept at different temperatures Tl and Tr and 
connected by a crossover domain of size £. Using the 
one-plate Green's function from Eq. (13) above, we can 
calculate the total stress on the plate as 



where VT = (1 — m)<5T /£, from which a Soret coefficient 
can be deduced as 



St 



A 



(23) 



Note that in this case the length scale £ appears only 
through the combination VT. 



C. discontinuous on the plate 

Finally, we consider a most general class of tempera- 
ture profiles as shown in Figs. 2(c) and (d), using the 
following example 



, T R + ST + e-"i 
— ' Tl + ST^ 



; z > 0, 
; z < 0, 



(24) 
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Due to the discontinuity of temperature at the bound- 
ary, the plate experiences the strongest form of force that 
reads 



n. 



k B 6T 



(25) 



to the leading order, where ST = (T R + ST+)- {T L + ST-) 
is the local temperature difference across the plate. Note 
that due to the singular form of the temperature profile, 
it is not possible to define a Soret coefficient in this case. 



V. TWO-PLATE GEOMETRY 

Here we consider the case of two parallel plates im- 
mersed in a nonuniform temperature profile. Wc com- 
mence by the equilibrium Casimir problem and then go 
ahead to generalize it to the nonequilibrium cases. We 
focus on two different types of behaviors corresponding to 
where the temperature profile is uniform or nonuniform 
in the exterior regions. 



T 











FIG. 3. Cross section of the temperature profile for the uni- 
form temperature case. The two plates are located at z = 
and z = H. 



The above stress has two terms; the first term is diver- 
gent and the second one is finite. The contribution from 
the left side reads 



IT 



lk B T f d 2 qduj I 



K 



(2tt) 3 Q++Q. 



(27) 



which exactly cancels the divergent part of Eq. (26), so 
that the total pressure on this plate is rendered finite as 



n, 



C(3) k B T 
8tt H 3 



(28) 



Note that the force is universal, and attractive (as it 
tends to push the plate at z — to the right). We can 
easily check that the other plate experiences the same 
amount of pressure in the reversed direction, and hence, 
the Newton's third law applies. As we will see below, 
this is a direct consequence of the fact that this force is 
mediated by equilibrium fluctuations. 




A. Equilibrium Casimir Force 

Here we derive the Fluctuation-induced force on the 
plates for the case when temperature profile is uniform 
everywhere, namely T(r, t) = To as shown in Fig. 3. The 
stress on the plate located at z = has two contributions 
coming from right and left sides. The contribution from 
the right is given as 



-k B T 



d 2 qdu! 



7 



(2tt) 3 
Q + [coth{Q+H) 



- 1] - Q_[coth(Q_iJ) 



2iu 



(26) 



FIG. 4. Cross section of the temperature profile for the 
non-uniform temperature case. The two plates are located at 
z — and z = H. (a) Temperature is uniform in the exte- 
rior, (b) Temperature is nonuniform in the exterior. In both 
cases the z-dependence of the temperature profile the interior 
is arbitrary. 

The cancellation of the divergent parts in Eqs. (26) 
and (27) is directly related to the assumption that tem- 
perature is continuous on the plates. If we considered 
the case where temperature is uniform but its value is, 
say, To in the region M and Tq in the region L, then in 
addition to the above Casimir force we would obtain a 
divergent part proportional to (To — Tq), similar to the 
single-plate forces discussed in Sec. IV C above. 
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B. Temperature profile independent of z in the 
exterior 



To extend the above calculations to the nonequilibrium 
situations, we first consider the case where temperature 
is independent of z in the exterior regions (R and L) but 
it is z-dcpcndent in the interior region (M), as shown in 
Fig. 4(a). Moreover, we also maintain the possibility of 
the dependence of temperature on time, as well as the 
in-plane coordinates r, and assume that the temperature 
is continuous on the plates to avoid the prescribed sin- 
gularities. 

We can calculate the two contributions to the local 
stress on the plate located at z — as 

f rfVrfa/ sinh(P + z') S inh(P_z') i(q . r _„ t) 
J (2tt) 3 sinh(P + i2>inh(P_i2") 

(29) 



and 



n = 



■yks f cPqduj 



I 



8K J (2,? 



1 



d 2 q'dcj' 



1 



(2tt) 3 P++P_ 



o*(q-r-wt) 



(30) 



where P± = \J \{q ± q') 2 + ^ UJ '~)> an ^ we have 

used the Fourier transform of temperature in time and in 
the in-plane coordinates. Due to the continuity of tem- 
perature, the divergent parts of the two contributions 
cancel each other as expected, and the net stress on the 
plate can be written as (see the Appendix) 



II(r,t) 



C(3) k B T(r,0,t) k B f d 2 qduj 



8tt 



H 3 



8irH 3 



T(q,0,w) 



(2tt)3 

e i(q-r-^) ) (3^ 



where 



f(u) = C(3) - -u 3/2 + 2u ln(l - e 2 ^ 7 ) 
+2V^ Li 2 (e 2 ^)-Li 3 (e 2 ^), 



(32) 



with Li„(z) = X^feLi z /k n being the polylogarithm func- 
tion. The function f(u) is analytic in the upper- half 
complex plane is asymptotes to £(3) everywhere on the 
semicircle at infinity. For small values of it, however, it 
has a power series expansion as: f(u) ~ u — |u 3 / 2 + 
£?i 2 + 0(u 5 / 2 ). An interesting feature of the above re- 
sult is that the stress on the plate is not sensitive to the 
^-dependence of the temperature profile in the interior 
region. 



In the case of uniform temperature Eq. (31) recovers 
the equilibrium Casimir force [Eq. (28)] because temper- 
ature will only have the zero frequency and wavevector 
component in Fourier space, and the correction term will 
vanish due to the factor of /(0) = 0. When temperature 
is modulated in time or in-plane space directions, the 
plate will experience local and instantaneous Casimir- 
like pressure as given by the first term in Eq. (31), while 
the second term yields corrections as 



n(r,t) 



C(3) fc B T(r,0,t) 



+ 



8tt 

k B 

16ttH V 2 



# 3 



+ K dt 



r(r,,,o,t), 



(33) 



to the leading order. Note that the above result is an ex- 
pansion in powers of the ratio -y-, where A is the typical 
wavelength set by the variations of the temperature in 
space or time, and is valid in the limit |<1. 

In the opposite limit of -y <c 1, the fluctuation- 
induced interactions are screened by the length scale A. 
This is reminiscent of the screening of the electric field for 
periodic charge distributions, where the screening length 
is also set by the period of the charge distribution. The 
screening manifests itself in the result of Eq. (31), where 
for large values of the wavevector or frequency the func- 
tion f(u) tends to C(3), thus leading to a cancellation of 
the first term. 

The above calculations can be repeated for the second 
plate to yield similar expressions for the stress, except 
that the value of the temperature will be replaced by its 
local values in the neighborhood of the second plate. In 
the general case when the two values are different, the 
two plates will experience different stresses and thus the 
third law of Newton will be violated. This effect is due 
to the nonequilibrium nature of the fluctuations, and has 
also been reported recently in the theory of depletion 
forces [9]. 

When temperature varies only with z, we can simplify 
the expressions for the stress on the left plate as 



Iltot(O) = 



C(3) k B T(0) 



8tt H 3 ' 
and correspondingly for the right plate as 

C(3) k B T(H) 



Tltot(H) = - 



8tt H 3 



(34) 



(35) 



Considering the fact that Newton's third law is violated, 
one can extract the pressure exerted on the system of the 
two plates as 



n 2 P C(3) fc B [T(ff)-T(0)] 
tot 8tt H 3 



(36) 



Note that the above result is reminiscent of Eq. (25), 
provided that we replace the microscopic length scale a 
by H , which is the effective thickness of the compound 
object. 



G 



C. Nonuniform temperature profile in the exterior 
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The fact that the fluctuation-induced interactions are 
independent of the interior temperature inhomogeneity in 
the z direction is a general result and does not depend on 
whether or not the exterior temperature is z-dependent. 
However, the form of the force that is acting on each 
plate depends on the neighboring exterior temperature 
profile. 

To study this class of temperature profiles, we con- 
sider the simple case where the temperature profile is 
independent of time and the in-plane space coordinates. 
We assume the following profile for the temperature 



T L + (T-T L ) e z /« 
T(z) = < arbitrary 

Tr + (T - T R ) e-(*" 



z < 0, 

< z < H, 

z > H, 



(37) 



as a typical example for the class of profiles depicted in 
Fig. 4(b). Note that we also assume that temperature 
is continuous on each plate. We can then calculate the 
stress on the left plate as 



n tot (o) 



C(3)fc B T , k B (T L -T) 



8tt m ' C3 ~ £a 2 
and correspondingly on the right plate as 



n to t(#) = - 



C(3)fcBT , . k B (T-T B 

8tt m 



c 3 - 



(38) 



(39) 



We can extract the pressure exerted on the system of the 
two plates as 



n 2 P 



<•?, 



k B {T L -T R ) 
£a 2 



(40) 



which is reminiscent of Eq. (25), and thus yields a Soret 
coefficient of the form given in Eq. (23) above. 



VI. CONCLUDING REMARKS 

In conclusion, we have shown that external objects 
located in a nonuniform temperature profile experi- 
ence nonuniversal fluctuation-induced forces that tend to 
move them towards the colder region, which is a manifes- 
tation of the so-called Soret effect. This effect is exam- 
ined in various temperature profiles, and it is shown that 
the behavior of the fluctuation induced forces strongly 
depends on: (I) whether the temperature is continuous 
across the object, (2) the asymptotic values of the tem- 
perature profile, and (3) temporal or spatial variations 
of temperature. In the case of two external objects im- 
mersed in a medium undergoing noncquilibrium thermal 
fluctuations, it is shown that the third law of Newton 
is violated in that the two objects experience different 
forces depending on where they are located. 
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APPENDIX A: INTEGRATION OVER A 
KERNEL 

In this Appendix, we outline some details of the cal- 
culations involved in Eqs. (29) and (30) above. We are 
dealing with the following integral: 



U(q,q!,w,z') 



dJ sinh(P + z')sinh(P_z') 
2tt sinh(P+tf)sinh(P_.H") : 



(Al) 



with P± being defined below Eq. (30). The integrand has 
two branch points, and their corresponding branch cuts, 
in the complex u>' plane. Now, by closing the integration 
contour around these cuts, we can directly check that the 
integral vanishes for z' < H, while it is divergent for the 
limiting case of z' = H . This behaviour is similar to a 
delta- function centered at z = H . We can then make the 
assignment: C/(q, q', oj, z') = u(q, q',w) S(z' — H) and 
try to find the function u. Note that if u is found to be a 
well-behaved finite function, our assignment will become 
rigorous. Integrating both sides of the above equation 
with respect to z' yields: 



«(q,q',w) = 



UK [ duj' P_ coth(P_#) 



7 
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2tt lj' - iK(q.q r )/j 



(A2) 



The integrand of Eq. (A2) has a pole and a branch cut 
as shown in Fig. 5. By closing the integration contour in 
the upper-half plane we obtain 



u(q,q',w) = -^y^ 2 + l' 2 ) +i ^ UJ 



x cothy^g 2 + q' 2 ) + i 7^ UJ , 
which completes the frequency integration. 



Complex CO Plane 



(A3) 




FIG. 5. Complex u/ plane and the contour that is used in 
the integration. The point P — iA"(q.q')/7 is the pole and 
BP — lo — I (q — q') is the branch point corresponding to 
the branch cut. 
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